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We show how Fermi liquid theory can be applied to ultra-cold Fermi gases, thereby expanding their "simula- 
tion" capabilities to a class of problems of interest to multiple physics sub-disciplines. We introduce procedures 
for measuring and calculating position dependent Landau parameters. This lays the ground work for addressing 
important controversial issues: (i) the suggestion that thermodynamically, the normal state of a unitary gas is 
indistinguishable from a Fermi liquid (ii) that a fermionic system with strong repulsive contact interactions is 
associated with either ferromagnetism or localization; this relates as well to 3 He and its p-wave superfluidity. 
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Normal Fermi liquid theory describes low temperature 
fermionic quantum matter which is adiabatically connected 
to its non-interacting counterpart |l|]. This theory, which 
has wide ranging implications for condensed and neutron star 
matter as well as nuclei and nuclear physics |2[], revolves 
around a quantification of Landau parameters. These are 
viewed as a collection of molecular field contributions char- 
acterizing the interaction between renormalized (quasi) par- 
ticles. Because Fermi liquid theory is a major tool for ad- 
dressing many-body systems and ultracold Fermi gases are a 
testing ground for counterpart theories, it is essential to es- 
tablish how to measure as well as calculate the associated 
Landau parameters. This leads to the goal of the present pa- 
per: to formulate the extension of Landau's important the- 
ory to an atomic trapped gas. Our work incorporates the lo- 
cal density approximation (LDA) which is used in both our 
experimentally-based and theoretical analyses. In the con- 
text of the latter, we present a methodology for calculating the 
Landau parameters, given a microscopic many body scheme. 

This paper addresses the relationship between Fermi liq- 
uid theory and two important situations: trapped gases in the 
presence of moderately strong attractive as well as repulsive 
contact interactions. For the former the ground state is a super- 
fluid in which the pairing strength can be continuously tuned 
from BCS to Bose Einstein condensation (BEC) yfl. For in- 
termediate strength attraction, associated with the "unitary" 
phase, of considerable interest is whether the normal state 
(above the transition temperature T c ) is a Fermi liquid or not. 
A body of evidence 14J-|6[] has arisen in support of the pre- 
diction [7] that somewhat above T c , an unpolarized, unitary 
gas contains pre-formed pairs associated with a normal state 
gap or pseudogap. If there is a gap in the fermionic excita- 
tion spectrum, these unitary gases would not qualify as Fermi 
liquids. However, this has become a more controversial point, 
recently when this observation was challenged 1 8] through fits 
to the measured temperature dependence of the pressure. The 
claim from Ref. [8] that "the low-temperature thermodynam- 
ics of the strongly interacting normal phase is well described 



by Fermi liquid theory" is both important and at odds with 
previous inferences 10-la 0] that Fermi liquid theory is inap- 
propriate at unitarity. 

By contrast, for the case of repulsive interactions, Fermi 
liquid theory should be widely applicable. Nevertheless, there 
are on-going controversies about the microscopic nature of 
such a Fermi liquid. In the cold gases this question is posed 
in the context of considering even stronger repulsion where 
the central issue is whether the ground state of the Hubbard 
gas model will be ferromagnetic Jiol [Till or localized lfl2ll . 
Here we argue that, in the Fermi liquid regime (at higher T, 
or weaker interactions so that there is no spontaneous symme- 
try breaking) one can use the behavior of the measured Fermi 
liquid parameters as a probe into the nature of the dominant 
incipient instability. Indeed, this debate in the cold gases is 
completely analogous to a related controversy in the proto- 
typical Fermi liquid of condensed matter physics: 3 He, which 
is thought to be similarly described by the repulsive Hubbard 
gas model. At issue is whether the Landau parameters are 
controlled by proximity to a ferromagnetic [13] or localiza- 
tion II 1 4TI instability. Further insight into 3 He should also be 
important to the cold gas community in additional contexts, 
since the p-wave superfluidity of 3 He at T = can serve as 
a template for addressing quite generally, non-s wave (e.g., 
d-wave) superfluidity which indirectly arises via repulsive on- 
site interactions. 

An understanding of precisely what does and what does 
not constitute a Fermi liquid is clearly lacking in the litera- 
ture of ultracold Fermi gases, where we have seen competing 
viewpoints expressed such as between Ref. |@] and Refs. Jj- 
0U>r between Refs. flolflTIl and Ref. 0, or Refs. and 
[14]. There has been some discussion of the basis of Fermi- 
liquid theory, albeit limited 11511 to the dynamical response in 
(less accessible) homogeneous systems. Fermi liquid theory 
is thought to be relevant to strongly polarized unitary gases, 
when superfluidity is destabilized. There are conjectures that 
a rapidly rotating unitary gas (such that superfluidity is sup- 
pressed) will also be a Fermi liquid [ 15], but this needs to be 
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systematically investigated since related analogous condensed 
matter systems show that pairing survives even when super- 
fluid coherence is destroyed fl. 

Because of these controversies and confusion, an impor- 
tant contribution of this paper is to establish the signatures of 
Fermi liquid and non-Fermi liquid behavior. In the process, 
we will also establish the thermodynamic signatures of pseu- 
dogap effects. Figure[T|is devoted to a comparison of a num- 
ber of Fermi liquid properties with those of a non-Fermi liq- 
uid. To represent the latter, we choose a unitary gas superfluid, 
for which there should be no disagreement that normal Fermi 
liquid theory fails. FigureQla) presents plots of the calculated 
density profiles for the Fermi liquid (FL) and non-Fermi liquid 
at the same temperature T = 0.5T C , where T c is the transition 
temperature of the unitary superfluid. For concreteness, in the 
FL case, we presume a repulsive contact interaction of mod- 
erate size and use Hartree-Fock theory (described below) to 
compute the profile. As an experimental calibration, for the 
unitary case both above and below T c , the profiles are in very 
good agreement [ 16] with the data. Despite their evident sim- 
ilarity, the FL and non Fermi liquid density profiles contain 
very different physics which must be carefully extracted. 

These differences are best accessed by studying derived 
physical properties. In this paper we focus on unpolarized 
Fermi gases. As a necessary first step, we compute the com- 
pressibility k, the spin susceptibility and the entropy from the 
density profiles. Within the LDA, the chemical potential is 
/i(r) = fi(r = 0) — (l/2)mw, 2 r r 2 , where m is the bare mass 
of the fermions and uj tr is the trap frequency. Then 
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which, in turn, allows one to measure k 



n~ 2 dn/dfi. 

From the profiles of a slightly polarized gas (with polariza- 
tion < 5%) one can arrive at the Pauli susceptibility Xs = 
dSn/dSfi ss Sn/Sfi. Here Sn — rt\ — 714. <C n and 
5(1 = — and one requires that /v and n a must be estab- 
lished separately. Measurements of \s have a firm basis in the 
past success of LDA approaches to polarized Fermi gases with 
attractive interactions [17]. One requires that the densities be 
small at the trap edge and that (for attractive interactions) low 
polarizations at T slightly away from zero are considered so 
that the phase separated state (with discontinuous densities) is 
not encountered. 

Techniques for extracting the entropy are directly based 
on those discussed elsewhere 111 811 . albeit modified somewhat. 
We consider two clouds at nearby temperatures T\ and T2 and 
label associated quantities of each cloud with subscripts 1 and 
2. Then by locating the points fii(Ti, r) = ^(Tz, r') on the 
two clouds fllSTI and using the Gibbs-Duhem relation, 
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Figure 1 : (Color online) Comparison of Fermi liquid (FL) and non- 
FL (unitary Fermi gas) thermodynamic behavior. In the FL case we 
consider a repulsive contact interaction with a = 0.2astoner- (a) 
Density profile of a unitary Fermi gas at T — 0.5T C (solid line) 
compared with that of the FL at the same T (dashed line). Here 
T c /Tf = 0.27. (b) Entropy profiles close to T c of unitary Fermi 
gases, slightly above (red) and slightly below (black) T c , and in the 
inset, at T = 0.5T C , as compared with that of the FL (dashed line) 
for the same T. (c) The spin susceptibility \s of a unitary Fermi gas 
(solid line, extracted from a 5% polarized gas with no phase sepa- 
ration) compared to that of the same Fermi liquid (dashed line) at 
T — 1.1T C . The inset shows counterpart plots for T = 0.5T C . Here 
fe_p(r = 0)astoner = it/2. Pseudogap effects in (b),(c) are reflected 
in the smooth evolution from above/below T c at unitarity. 



We note that the extracted s(r) agrees very well with the ex- 
act s(r) provided one uses the averaged T. With care, one 
should be able to choose the separation between T\ and T2 
large enough so that noise effects do not affect the extraction 
of s(r). 

We show in the inset to Fig. Lb) a plot of s(r) calcu- 
lated d 1 911 at unitarity as compared with the FL case, for 
T/T c — 0.5. In the main body of Fig. |TJb) we compare 
this with the counterpart unitary plots just below the transi- 
tion temperature (T/T c = 0.9) and in the normal phase at 
T/T c = 1.1. The total entropy involves a bosonic as well 
as a fermionic contribution 11911 . A salient feature of s(r) in 
these unitary superfluids is that it is spatially non-mono tonic. 
This behavior is associated with a pairing gap which is largest 
at the trap center and decreasing with increasing r. That this 
non-monotonicity in s(r) is found at T/T c = 0.9 suggests 
the pairing gap is still large at these high temperatures; be- 
cause the transition is second order this pairing gap must per- 
sist into the normal phase, as a "pseudogap" O. I9I I20I1 asso- 
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ciated with non-condensed pairs. Indeed, one sees that the 
curves for T/T c = 1.1 exhibit an (albeit, slightly) anomalous 
non-mono tonic entropy profile. 

Similar non-monotonic behavior can also be found in the 
spin susceptibility \s of an unpolarized unitary Fermi gas 
which is plotted in Fig.QJc) for T/T c = 1.1 (and in the inset 
at T/T c = 0.5) compared to that of the same Fermi liquid as 
discussed in Fig.[Tfa). These two non-monotonicities contrast 
with their counterparts for the Fermi liquid case shown in the 
figures. The suppression of \s at the trap center of the normal 
unitary gas can similarly be traced to singlet non-condensed 
pairs 01711 . or pseudogap effects. 

The behavior of the normal unitary gas, thus, shows a con- 
tinuous evolution from the superfluid properties ( say, mea- 
sured at T/T c rs 0.5) to above T c . This is the sense in which 
pseudogap effects (manifestly inconsistent with Fermi liquid 
theory) will be evident in thermodynamics. This is in contrast 
to the claim, based on a quadratic temperature dependence of 
the pressure [8] that "the normal phase of the unitary Fermi 
gas., [has] thermodynamic properties [which are] .. well de- 
scribed by Fermi liquid theory, unlike high-T c copper oxides." 
The same quadratic signature [8] can also be found theoreti- 
cally for a range of T in strict BCS theory, where FL theory 
obviously fails. Thus, we infer that power law dependences in 
thermodynamical variables are not valid indicators of the ap- 
plicability of Fermi liquid theory. Rather, non-monotonicity 
in s(r) and Xs profiles, which reflects the T < T c physics, 
is the more relevant feature, a consequence but not a proof of 
the existence of a pseudogap [0-01 and the concomitant failure 
of Fermi liquid theory. Interestingly, a number of the theories 
with which these experiments [8] are favorably compared also 
contain pseudogap effects. 

In proper Fermi liquids the compressibility, spin suscep- 
tibility and entropy enable one to experimentally extract the 
spatially dependent Landau parameters [ 1] and effective mass 
to*. In condensed matter Fermi liquids ra* is most readily de- 
duced from the linear (in temperature) term in the specific heat 
U]. In cold gases, we argue that it is best obtained from the 
linear in temperature term in the local entropy density s(r): 



s(r) = 



ra*kp(r)T 



(3) 



Here we set h = 1 and ks = 1. Fermi liquid theory and 
Eq. (O are precise only if T <C Tp. As the local T F (r) de- 
creases, this procedure will fail towards the trap edge. Thus 
expressions including only leading order T corrections in s(r) 
(or equivalently T 2 terms in P or \x, etc.) must be constrained 
to T <C Tp(f). Moreover, in a Galilean invariant system U, 
Ff satisfies 
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(4) 



and these arguments apply to a trapped gas within the LDA. 
The Landau parameter Fq can be obtained from 



dn ( dno \ m* /ra 
dp ~ \~dp) 1 + Ffi ' 



Here da /d[j, = - 2k( 9 /( w )/ 9u; )U=ft 2 fe 2 /2m-Ai is the 
known corresponding quantity for a non-interacting Fermi gas 
at the same chemical potential and f(x) is the Fermi distribu- 
tion function. Given dn/dn{r) and m*(r)/m, Fq (r) is then 
known. Finally, the third Landau parameter to be quantified is 
Fq, which is related to the Pauli susceptibility 
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Here XsO = dSno/dSfi is the spin susceptibility of the non- 
interacting Fermi gas. 

We consider now a repulsive contact interaction in order 
to gain more understanding of the stable phases of the Hub- 
bard Fermi gas Hamiltonian. This bears on recent 
and longstanding controversies lfl3l [l4ll in the literature. If 
the repulsion is not too large, one expects the system to be 
in a Fermi liquid phase. Our methodology for computing 
the associated Landau parameters is based on first comput- 
ing the trap profiles n a (r) from a microscopic theory and sec- 
ondly following the experimental protocols outlined above to 
determine the LDA-derived Ff, Fq, and Fq. The specifics 
of a microscopic theory are tested by comparing the calcu- 
lated Landau parameters with their experimental counterparts. 
In LDA the density is n a (r) = J2w,u n G (T (S cr (k, ui n> n(r)). 
Here G a is the Green's function, u n is the fermion Matsub- 
ara frequency, and the fermion self-energy Y, a (k, w„, //(?")) is 
chosen to correspond, for example, to a nearly ferromagnetic 
IU3I1 or nearly localized d 1 411 theory. 

While it is relatively straightforward to use more sophisti- 
cated approaches [13], in order to focus on the central physics, 
we apply Hartree-Fock theory 111 ll |2lll . Here Y, a = gn s , 
where g = Ama^a/m, a = —a, and a is the two-body s-wave 
scattering length. The trap profiles correspond to n a (r) — 

/ (t^T + 9 n s{r) - n a {r)\ . We consider a/a sto „ er = 
0.6 and 0.9 at T/T F = 0.02 and 0.01, where a S toner is the 
critical scattering length (Stoner instability) at which Xs di- 
verges. From the profiles, one obtains m*, (dn/dfi), and Xs- 
The corresponding spatial dependences associated with each 
of these are shown on the left column of Fig. [2] where we fo- 
cus on the results at the trap center. Using Eqs. ©, and 
the three Landau parameters F* , Fq, and Fq are obtained 
and plotted in the right column of Fig. [2] We see that m* /m 
is only slightly larger than 1 so F( is small. For repulsive in- 
teractions, Fq is expected to be positive and Fq is expected to 
be negative, as is consistent with Fig.|2]panels (e) and (f). 

Interestingly, had one naively applied these same FL pro- 
cedures to the attractive case (where we argue because of the 
pseudogap, Fermi liquid theory is not appropriate) one would 
find that Fq is positive and Fq is negative. Therefore, even 
though the trap profiles for the particle density are not qualita- 
tively different for both cases, the derived Landau parameters 
make it possible to distinguish between very different physical 
situations. 

One can see from Fig. |2|d) and (f) that as a — )• astoner 
Fq — >• — 1 while Ff and the effective mass remain finite. In 
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pseudogap effects is comparable to that in Refs. |H , 22 ] . In this 
paper, following Ref. 10 we have also mapped out a new direc- 
tion for the field of ultra-cold trapped Fermi gases which, via 
simulations of the (repulsive) Fermi Hubbard gas, can serve to 
establish where ferromagnetism, localization and p wave su- 
perfluidity are stable. This should be facilitated by analyzing 
simultaneous measurements and calculations of appropriate 
Landau parameters. 
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DMR-02 13745. We thank Q. J. Chen for providing thermo- 
dynamical plots. C.C.C. acknowledges the support of the U.S. 
Department of Energy through the LANL/LDRD Program. 



0.5 



Figure 2: (Color online) The extracted (a) m*/m, (b) 
(dn/dfi) I (dno /dp), and (c) Xa/XsO and the Landau parameters (d) 
F", (e) F§, and (f) Fg . The input data are from the Hartree-Fock 
theory with a/astoner = 0.6 (black dashed line) and 0.9 (red solid 
line) at T/Tf = 0.02 (whose n(r) and p(r) are shown as insets) 
and 0.01. 



contrast, for an almost localized liquid [12], F{ or m*/m will 
diverge at the counterpart critical interaction strength. While 
this leads to a divergence in \s, F§ does not approach —1 
lfl4ll . For an unknown ground state, this Fermi liquid analysis 
in the limit that the interaction strength is slightly reduced, 
can be used to provide important hints about the nature of the 
more strongly correlated phase. 

Thus, future measurements of the Landau parameters as- 
sociated with the repulsive Hubbard gas studied in Ref. llOl 
should be able to shed light on a current controversy 111211 
over whether a ferromagnetic ground state has been observed. 
There are, however, some complications in these experiments. 
The Feshbach resonance which tunes the repulsive interac- 
tions between fermions is associated with a molecular bound 
state lying below the particle continuum which is populated 
ifioll at the 25% level. These molecules derive from molecule- 
particle or molecule-molecule scattering which may lead to 
systematic errors in determining dn/dp and \s- To resolve 
this issue, one has to either reduce the population of molecules 
or develop a more sophisticated theory to take them into ac- 
count. 

An important contribution of this paper has been to estab- 
lish how, despite previous claims [81], pseudogap effects do 
enter into thermodynamical measurements. When a pseudo- 
gap is present, thermodynamical variables above T c will show 
a smooth evolution from their behavior deep in the superfluid 
phase; this contrasts with BCS theory, where the normal state 
thermodynamics do not reflect their superfluid counterparts. 
While there will be a thermodynamical feature of the phase 
transition ||9|] one should see this smooth evolution in anoma- 
lous spatial dependences for example of the spin susceptibil- 
ity. The experimental complexity needed to search for these 
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